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ABSTRACT 

We argue that the X-ray and UV flux illuminating the parsec-scale accretion disk around luminous 
active galactic nuclei (AGN) is super-Eddington with respect to the local far-infrared dust opacity. 
The far infrared opacity may be larger than in the interstellar medium of the Milky Way due to 
a combination of supersolar metallicity and the growth of dust grains in the dense accretion disk. 
Because of the irradiating flux, the outer accretion disk puffs up with a vertical thickness h ~ R. This 
provides a mechanism for generating a geometrically thick obscuring region from an intrinsically thin 
disk. We find obscuring columns 
Subject headings: galaxies: active - galaxies: Seyfert 
dust, extinction 



10 22 — 10 23 cm 2 , in reasonable agreement with observations. 



galaxies: nuclei - accretion, accretion disks 
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1. INTRODUCTION 

The Seyfert unification model postulates the existence 
of a dusty geometrically thick obscuring region (the 
"torus" ) in order to account for many of the observed dif- 
ferences among various classes of AGN (Antonucci 1993). 
The large infrared fluxes from quasars also provide evi- 
dence for obscuration of the central source, with covering 
fractions of 10-50% (Sanders et al. 1989). Recent obser- 
vations with SDSS have directly detected a large number 
of Type 2 (obscured) quasars (e.g. Zakamska et al. 2003, 
2005; Ptak et al. 2006). 

The origin of the geometrically thick obscuring mate- 
rial in AGN remains uncertain. The host galaxy may 
provide the obscuration as in the case of MCG-6-30-15 
(Ballantyne, Weingartner, & Murray 2003). A warped 
disk is directly implicated in the low-luminosity AGN 
NGC 4258 (Greenhill et al. 1995) and may be relevant 
more generally (Sanders et al. 1989). Alternatively, the 
obscuring material may arise in an outflow from the un- 
derlying accretion disk (e.g., Konigl & Kartje 1994). En- 
ergy injection by star formation (e.g., Wada & Norman 
2002; Thompson et al. 2005) or viscous stirring in a 
clumpy accretion disk (Krolik & Begelman 1988) may 
also help generate the large random motions necessary 
to maintain h ~ r. Lastly, Pier & Krolik (1992) pointed 
out that, for luminous AGN, radiation pressure from the 
central AGN could help maintain the vertical thickness 
of obscuring material. 

In this paper, we present a model for how geometrically 
thick obscuring material might arise from X-ray illumi- 
nation incident on an initially thin AGN disk (§[2|). This 
model is complementary to that of Pier & Krolik (1992); 
we focus on how X-ray irradiation can thicken an under- 
lying thin disk, while they highlighted the importance 
of radiation pressure support for already geometrically 
thick material. In § |3j we argue that the far infrared 
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opacity in AGN disks is larger than that in the interstel- 
lar medium (ISM) of the Milky Way, which increases the 
dynamical importance of X-ray and UV irradiation. The 
dust opacity is enhanced due to 1. supersolar metallicity 
and 2. larger dust grains. We discuss the evolution of 
dust grains in AGN disks via accretion of metals from 
the gas and dust grain coagulation. We show that the 
conditions in a thin AGN disk appear conducive to grain 
growth (as also suggested by, e.g., Laor & Draine 1993 
and Maiolino et al. 2001ab). Finally, in § H we discuss 
our results and their observational implications. 



2. X-RAY HEATED AGN DISKS 

Here we consider the impact of X-ray illumination on 
the structure of AGN disks. Our physical picture is as 
follows (see Fig. [IJ. UV and X-ray radiation from the 
central source are incident at a grazing angle on a region 
of a thin disk which has a thickness h/R. Due to the high 
opacity of dust to UV photons, the incident UV flux is 
absorbed well above the IR photosphere of the disk. UV 
irradiation thus tends to make the disk more isothermal 
and does not provide an internal radiation pressure gra- 
dient which can help support the disk. X-rays in the 10- 
100 keV range, on the other hand, are electron-scattered 
in the outer atmosphere of the disk with approximately 
half of the X-rays being downscattered. If the underly- 
ing accretion disk is sufficiently Compton thick, most of 
the downscattered X-rays lose their energy via Comp- 
ton scattering below the IR photosphere. This energy 
is then re-radiated in the IR, diffuses towards the sur- 
face of the disk, and provides a radiation pressure gra- 
dient which is dynamically important if the energy de- 
position via X-rays is sufficiently high. We show that as 
a result of this X-ray irradiation, the entire disk cannot 
remain thin. Instead, the outer atmosphere puffs up to 
h ~ R producing significant obscuration along most lines 
of sight. As discussed later in §4, our calculation does 
not quantitatively describe the final structure of the re- 
sulting geometrically thick atmosphere, which requires a 
multi-dimensional calculation that incorporates both UV 
and X-ray irradiation. 

We now calculate the critical luminosity so that X-ray 
irradiation becomes dynamically important. The inci- 
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X-ray and UV radiation 




Fig. 1. — Schematic Diagram of an X-ray Illuminated AGN Disk 



dent X-ray flux on the initially thin disk is 
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where Lx is the X-ray luminosity from the central source 
and R is the radius of the disk. The X-ray opacity is 
provided by electron scattering and photoionization. For 
X-ray energies above w 10 keV, electron scattering dom- 
inates. At these energies, an x-ray photon incident on 
a disk with scale height h/R is scattered at a column 
of 1/ K e (h/R), where K e is the electron scattering opac- 
ity. Since the scattering is roughly isotropic, half of the 
X-rays are scattered out to infinity and half are scat- 
tered down into the disk. The half that are scattered 
downward heat the disk via photoionization and Comp- 
ton scattering. We adopt the photoionization cross sec- 
tion from Maloney, Hollcnbach, & Tielens (1996) which 

gives cr = 4.4xl0~ 22 (E/l keV)~ 8/3 cm" 2 for E > 7 keV, 
where E is the energy of the X-ray photon. We calculate 
the efficiency of this heating for an incident spectrum 
typical of AGN (photon index a — 1.7) as a function 
of depth and show the results in Figure O The total 
energy absorbed by the disk is a function of its Comp- 
ton thickness, but above a total column, y tot = J pdz, 
of 20gcm~ 2 , the heating saturates. For these columns, 
w 2/3 of the energy of the downward scattered X-rays 
heats the disk, with most of the energy being deposited in 
the first few Thomson depths. Models of accretion disks 
in luminous AGN predict surface densities of ~ 10 — 10 3 g 
cm~ 2 on parsec scales (Sirko & Goodman 2003; Thomp- 
son et al. 2005), so X-ray heating is expected to be quite 
important. These X-rays are reprocessed into the in- 
frared where the opacity is substantial. 

Under the conditions of interest, the radiation pressure 
gradient typically exceeds the gas pressure gradient near 
the IR photosphere. In order to support the photosphere 
to a height h, the IR flux, Fir, through the disk must 
satisfy the dust Eddington limit (Fir w -FEdd,d) or 
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where Kd is the dust IR opacity. Since the net flux radi- 
ated by the IR photosphere is due to the fraction of the 
incident X-ray flux deposited beneath the photosphere, 
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Fig. 2. — Local heating rate (H x /p) per unit mass and inte- 
grated heating rate (j(H x /p)dy) as a function of column y (mea- 
sured downward into the disk) for a slab irradiated by a power-law 
X-ray spectrum with a photon index of 1.7 and a high energy cutoff 
of 100 keV. Arrows indicate the vertical axis associated with each 
solid line. The heating rates are normalized to the total x-ray flux 
incident on the column (F x ) so that (J H x dy) / F x is the fraction of 
the incident flux that has been absorbed at a given depth y. The 
total heating rate, Compton scattering heating, and photoioniza- 
tion heating are given by the solid-lines, dotted lines and dashed 
lines, respectively. 



we also have 
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where e w 1/3 is the efficiency with which X-ray flux is 
absorbed deep in the disk, which we estimate from the 
half of the X-rays that are downward scattered and the 
two-thirds of this energy that is absorbed. For a warped 
disk, all of the incident X-ray may heat the disk, not just 
the scattered radiation. Equating equations ([2]) and (|3|), 
we find a critical luminosity 



_ 1 4irGM B nc _ _ x K e 

^X, crit — £ — e -^Edd , 



(4) 



where L-Edd = 4,ttGMbh_c/ K e is the electron scattering 
Eddington limit and Mbh is the mass of the central black 
hole. For Lx > L-x.crit, the disk is locally unstable to 
puffing up because the energy deposited below the IR 
photosphere exceeds the local dust Eddington limit. 

The critical luminosity Lx, crit is sensitive to the details 
of the dust opacity in AGN disks through n e /nd- For the 
conditions of interest T ~ 300 — 1000, K and normal ISM 
dust is estimated to have a Rosseland mean opacity be- 
tween Kd ~ 5 (Semenov et al. 2003) and Kd ~ 15 cm 2 g _1 
(Drainc 2003). In the next section, we argue that the 
far infrared opacity in the near- AGN environment may 
be as large as Kd ~ 50cm 2 g _1 due to increased metal- 
licity and grain growth in AGN disks. In this case, 
Kef^d ~ 0.01 which in turn implies ix,crit/-^Edd ~ 0.03. 
Bright Seyferts and quasars have X-ray luminosities that 
are roughly 0.01-0.1 of the (electron-scattering) Edding- 
ton luminosity, so a large number of systems are near the 
critical luminosity needed to have X-ray heating signif- 
icantly affect the structure of a thin dusty parsec-scale 
disk (§4). 
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To understand the impact of X-ray irradiation in more 
detail, we calculate the vertical structure of an irradiated 
thin disk. We begin with the equation of hydrostatic 
balance, 



dP t 



dz 



= -p\9 



K d F 



and radiative balance, 
dP, 



rad 



KdFp 



dz 



(5) 
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where P gas = pkBT/m p is the gas pressure, T is the 
temperature, P ra d = aT 4 /3 is the radiation pressure, 
F is the flux, z is the height defined from the mid- 
plane, g — 2ttGYi + il 2 z is the local gravitational ac- 
celeration, and = */ GMbh/R 3 is the local Keple- 
rian frequency. Note we define the disk surface den- 
sity £ = J pdz = at the midplane here and is de- 
fined upward. For the Rosseland mean opacity of dust, 
we assume a simple model for numerical convenience: 
we take a constant value of Kd = 50 cm 2 g -1 between 
T = 100 K and T = T sub = 1500 K (see eq. Q7]). We 
model the rapid drop in opacity as the dust sublimates 
using K d oc exp(— (T- T suh )/AT) for T > T sub , where 
AT w 100 K defines the rough width of the transition. 
None of the qualitative features of our model depend on 
the choice of Kd — 50 cm 2 g" 1 at low T, although the 
quantitative importance of irradiation does depend on 
the absolute normalization of Kd (as discussed above in 
the context of Lx,crit)- 

The impact of X-ray heating on the structure of the 
disk can be modeled as a source term in the equation for 
the flux: 

dF 

=- = H{z) (7) 
dz 

where H(z) is the heating function shown in Figure[2] To 
calculate the structure of the irradiated disk, we integrate 
equations (JSJ), ©, and ((7J from z = to z = z , where 
zq is the IR photosphere. We enforce F(z — 0) = at 
the midplane and the photospheric condition at z = zq, 
t sa Kphp — 2/3, where r is the optical depth taken 
from infinity and h p is the local scale height of the pho- 
tosphere. Numerically we specify the height of the IR 
photosphere and compute the required incident X-ray 
flux. 

Figure [3] shows the scale height h/R of an irradiated 
disk as a function of the X-ray luminosity of the central 
source. The X-ray luminosity is normalized to Lx, exit 
from equation (0J . Two representative vertical structure 
models (h/R = 0.01 and 0.5) are shown in Figure [U In 
our calculations, we determine the downscattered X-ray 
flux, Fx,ds, needed to support a disk with scale height 
h/R. In Figured this is interpreted as an X-ray luminos- 
ity of the central source using Lx = 2x4irR 2 Fx,ds(R/h). 
The factor of 2 is due to the fraction of the X-rays down- 
scattered in the atmosphere of the disk and the factor of 
R/h is due to the grazing incidence of X-rays (eq. [1]). 
The results in Figures [3] and H] are for a disk with a surface 
density of 50 gem -2 around a 10 7 M Q BH at R = 0.2 
pc. The total column in this case is such that nearly all 
of the downscattered X-rays are absorbed (c.f. Figure 
[2|). The structure of the surface layers of the disk (e.g., 
the photospheric height h) are independent of the total 
column once it exceeds « 20 — 30 g cm -2 . Note that 
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Fig. 3. — Scale height h/R of an X-ray irradiated disk as a func- 
tion of the X-ray luminosity of the central source, for a disk with a 
surface density of 50gcm~ - 2 around a 10 7 M e BH at R = 0.2 pc. 
The X-ray luminosity is normalized to Lx.crit from equation Q. 
We calculate vertical structure models and determine the down- 
scattered X-ray flux Fx ds needed to support a disk with a scale 
height h/R. The required X-ray luminosity of the central source 
is then estimated to be Lx = 2 X 4irR 2 Fx,ds{R/ 'h) (see the text 
for details). Note how a factor of 2-3 increase in Lx increases the 
scale-height of the disk by a factor of ~ 100, from h/R w 0.01 to 
h/R « 1. 

for the solutions shown in Figures [3] and there are 
no internal sources of heat, such as would be provided 
by viscous stresses or star formation; including such in- 
ternal sources of heat would lead a minimum h/R even 
in the absence of irradiation (unlike in our present solu- 
tions in which h/R — > as Lx — > 0). This would not 
significantly change the properties of our solutions when 

L X ~ £*X,crit- 

Figure [3] shows that for Lx <C ix.crit, the irradiating 
X-rays have little effect on the structure of the disk, even 
in the photospheric layers: the disk remains cold and thin 
and is gas pressure dominated throughout. X-ray irradi- 
ation can, however, still affect the chemistry of the disk 
in this limit, and can lead to conditions conducive to wa- 
ter masing ala NGC 4258 (Maloney et al. 1996). For low 
accretion rates, the surface density of an AGN disk can 
also be sufficiently small that X-ray irradiation leads to 
a transition from molecular to atomic gas (e.g., Neufcld, 
Maloney & Conger 1994). The gas surface densities in 
bright AGN are, however, too large for this to occur. 

For Lx — > Lx, crit; Figure [3] shows that the disk puffs 
up to h ~ R. Even in this limit, however, most of the 
mass remains relatively unaffected by the X-ray irradi- 
ation (e.g., the half-mass heights for the models in the 
left and right panels of Fig. H are 0.002 R and 0.02 
R, respectively). The transition between X-ray irradi- 
ation having relatively little dynamical effect to forcing 
the photosphere up to h ~ R is quite abrupt, with h/R 
increasing by a factor of ~ 100 over only a factor of 2.5 in 
Lx ■ Physically, this is because at the temperatures of in- 
terest a gas pressure dominated disk has h <C R and only 
when Lx ~ Lx, crit does radiation pressure contribute 
significantly to the pressure support. Observationally, 
the rapid change in h/R in Figure [3] would correspond to 
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Fig. 4. — Vertical structure of an X-ray irradiated disk around 
a 10 7 M o BH at R = 0.2 pc. The incident X-ray flux is ad- 
justed to give h/R = 0.01 (left) and 0.5 (right). In the upper 
panel we plot the temperature as a function of height. In the 
lower panel we plot the density (solid line), flux (dashed line) and 
column (dotted line) normalized to the midplanc density (p m id), 
photospheric Eddington flux (Fgdd d = GMh/{R i n^)), and to- 
tal column Stot = 50gcm~ 2 , respectively. The column, £, is 
defined to be zero at the midplane. For h/R = 0.01, the disk 
is gas pressure dominated. The h/R = 0.5 solution, on the 
other hand, is radiation pressure dominated (F/Fgdd d ~ 1) above 
z/h ~ 0.05. The half-mass heights (midplane densities) are 0.002 
R (n ss 1.5 X 10 10 cm- 3 ) and 0.02 R (n K 2 X 10 9 cm" 3 ) for 
h/R = 0.01 and h/R = 0.5, respectively. The incident X-ray flux 
is F x ~ 3.3 X 10 6 ergs cm" 2 s" 1 (L x /L Edd « 0.01) for h/R = 0.01 
and F x ~ 7 x 10 6 ergs cm" 2 s" 1 (L x /L Edd « 0.03) for h/R = 0.5. 

a dramatic change in the nuclear obscuration on parsec 
scales for systems with Lx > Lx.crit- 

For an irradiating flux of Lx > L x , C rit, the atmosphere 
of the irradiated disk becomes geometrically thick and 
the outer layers will be blown away by a combination of 
the X-ray and UV radiation. Though a full wind calcu- 
lation is needed to understand the dynamics of such a 
flow in detail, we can roughly estimate the column den- 
sity through the unbound material. Since the wind is 
momentum driven by a super-Eddington flux, the maxi- 
mum mass outflow rate is given by momentum balance, 
M max t> csc = L/c, where v csc is the escape velocity from 
radius R in the disk, which gives 
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10 44 ergs s- 



(8) 

where R x = (i?/lpc) and M 7 = (M BH /10 7 M Q ). The 
obscuring column through such a flow is 
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For typical values of L/L^dd in luminous AGN, the col- 
umn through the outflowing wind is Compton-thin with 
N H « 10 22 - 10 23 cnT 2 . Equations (J5J) and © assume 
that all of the luminosity of the central source contributes 
to driving an outflow. This is not correct for the X- 
ray flux, since the resulting column is Compton thin. 
However, once the photosphere of the disk is moderately 
inflated by X-ray irradiation, the UV flux will become 
dynamically important and will also contribute to un- 
binding the surface layers of the disk. Thus the above 
estimates of the mass loss rate and column are likely to 
be reasonable. 

3. DUST GRAIN GROWTH IN AGN DISKS 

As discussed in the previous section, the impact of 
the X-ray and UV radiation from a central AGN on 
its surrounding accretion disk depends sensitively on 



the FIR opacity of the disk. For temperatures of ~ 
300 — 10 3 K, typical models of ISM dust have opaci- 
ties of 5 — 15 cm g . We now argue that this opacity 
may be enhanced in the near- AGN environment due to 
1. supersolar metallicity and 2. larger dust grains. 

If most metals are locked up into dust, enhancements 
in the metallicity increase the opacity as Z/Zq. From 
quasar observations, the broad line region appears to 
have metallicities that are a few to ten times solar 
(Hamann & Ferland 1999; Baldwin et al. 2003; Diet- 
rich et al. 2003; Nagao, Marconi, & Maiolino 2006). 
Note that these metallicities probe just the central nu- 
clear regions of AGN, whose chemical evolution might 
be significantly different from the host galaxy (Hamann 
& Ferland 1999). In addition recent work on the narrow 
line region of Seyferts and AGN outflows show supersolar 
metallicity (Groves, Heckman & Kauffmann 2006; Arav 
et. al. 2006). For Draine's (2003) Rosseland mean opac- 
ity of Kd w 15 cm 2 g -1 , a metallicity enhancement at the 
level of 2-4 would increase Kd to « 30 — 60 cm 2 g _1 . This 
corresponds to n e / Kd ~ 10~ 2 , in which case AGN with 
X-ray luminosities > 0.01 — 0.03 LEdd wm have an X-ray 
flux that is super-Eddington with respect to the dust on 
parsec scales. 

An independent argument for enhanced FIR opacity in 
small-scale AGN disks is that the dust grains may pref- 
erentially be larger than in the local ISM. We begin by 
considering the observational lines of evidence for (and 
against) this conclusion. Recent detections of the 10 fim 
silicate emission feature in AGN find that the emission 
feature is broadened and shifted to slightly longer wave- 
lengths than in the ISM (Siebenmorgen et. al. 2005; Hao 
et. al. 2005; Sturm et al. 2005). This is analogous to the 
emission features in Herbig Ae/Be systems (Bouwman 
et al 2001) and may indicate different dust properties, in 
particular larger or more amorphous grains. On the other 
hand, the silicate absorption features in Type 2 AGN do 
not appear to require unusual dust properties, but may 
require a patchy distribution of dust (Roche et al. 1991; 
Hao et al. 2007; Levenson et al. 2007). Additional con- 
traints on the properties of dust in the vicinity of AGN 
can be provided by extinction measurements. Maiolino 
et al. (2001ab) used such measurements to argue that 
dust grains near AGN may be significantly larger than in- 
terstellar medium grains. Similarly, Gaskel et al. (2004) 
argue that the reddening law near AGN is substantially 
different from any standard reddening law and that the 
minimum size for dust grains in the near- AGN environ- 
ment is a ~ 0.3/Ltm due to a lack of reddening in the UV. 
Using SDSS data on a large number of AGN, however, 
Hopkins et al. (2004) argue against such a Gaskel-like 
reddening law, and in favor of SMC-like dust. 

Although the observational arguments for larger dust 
grains near AGN are by no means conclusive, there are 
theoretical arguments for this conclusion as well. To 
show this, we consider the evolution of dust grains in 
AGN disks via accretion of metals from the gas and dust 
grain coagulation. Lacking a detailed model for the ISM 
on parsec scales near AGN, we make two simplifying as- 
sumptions. First, we assume that grain growth is im- 
portant if it occurs on a timescale less than the local 

1 3 /2 — 1 /2 

dynamical time td yn ~ ^ ~ 5000-R/ M 7 ' yrs. The 
dynamical time is comparable to the eddy turnover time 
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on the outer-scale, which is ~ the timescale on which 
dust would be destroyed by shocks if turbulence is super- 
sonic in the disk. Our second simplifying assumption is 
that we estimate the gas density in the disk p g in terms 
of the Toomre Q parameter for gravitational stability 



via p g = M bh Q- 1 /(2ttR 3 



6 x 10 7 Q- 1 M r Rl 



i 



cm . We assume Q ~ 1, as suggested by theoretical 
arguments (e.g., Sirko & Goodman 2003; Thompson et 
al. 2005) and as is observed in a diverse range of galax- 
ies from local spirals (e.g., Martin & Kennicutt 2001) 
to luminous starbursts (e.g., Downes & Solomon 1998). 
Although the ISM in the vicinity of AGN undoubtedly 
consists of multiple phases, Q ~ 1 should provide a rea- 
sonable guide to the mean ISM densities. The gas den- 
sities implied by Q ~ 1 are also consistent with the den- 
sities of ps 10 7 — 10 10 cm -3 required to account for the 
observed water maser emission from the central parsecs 
of a number of AGN (see, e.g., Lo 2005 for a review). 

We begin by considering the timescale for a dust grain 
to grow appreciably in mass via accretion of gas-phase 
metals, which is 



4 apd 
3 v m p m 



, nn aiPd.zRlQ „. 
100 — yrs < i dyn 



T 3 1/2 M 7 



1/2 



(10) 

where m d = ( o ci 47ra 3 /3 is the mass of the dust grain, a 

na 2 is the dust-gas collision cross section, v m « 1 T, 
km s _1 is the thermal velocity of the metals at temper- 
ature 10 3 T 3 K (assuming A ~ 20), pd — 3pd,3gcm~ 3 is 
the dust grain density, a — aipm is the radius of the 
dust grain, and p m ps 0.01p g is the mass density in met- 
als. Note that we take a standard dust grain density of 
3gcm~ 3 . Equation (flQ| shows that dust grains can in- 
crease their mass on a timescale significantly less than 
the local dynamical time. Thus even if some grains (or 
their mantles) are destroyed, the free metals accrete onto 
remaining dust quite rapidly. 

The high gas densities in AGN disks also make them 
a promising environment for grain coagulation, as is in- 
ferred to occur in protostellar disks (van Boekel et al. 
2003). To show this, we estimate the characteristic ve- 
locity dispersion of dust grains and show that a signifi- 
cant fraction of small dust grains stick on impact. The 
timescale for coagulation is of order the dynamical time 
so this process occurs quickly. 

We begin by estimating the velocity dispersion of 
grains in a turbulent medium following the argument 
given by Draine (1985). The velocity dispersion of grains 
is determined by the scale on which the eddy turnover 
time, teddy, is comparable to the aerodynamic drag 
timescale, idrag, for grains in gas. The latter timescale is 
given by 



tdi 



m d 

Pg(TC s 



8 x 10 6 pd,3-RiMf ^3 1/2 ai Q s, (11) 



where c s is the thermal velocity of the gas. On the largest 
scales, the accretion disk possesses a turbulent velocity 
approaching h/R times the orbital velocity, u or b- As in 
the ISM, these velocities are likely supersonic compared 
to the gas sound speed and so we assume a supersonic 

1/2 

cascade with v e ddy ~ (h/R)v OI b (l c ddy/h) , where l e ddy 
is the size scale of the eddy. The eddy turnover time, 



teddy, is given by 



^eddy ^ ^dyn 



^eddy 



1/2 



(12) 



As turbulence cascades to smaller and smaller scales, 
the eddy velocities decrease. Eventually, f ddy ~ c s 
and turbulence transitions from supersonic to subsonic. 
We define the eddy turnover time at this scale (when 
Veddy = c s ) as the sonic time, t S onic, which is given by 



tdy 



Res 

' h for 



100T 3 1/2 Aff ^ 



yrs (13) 



The sonic time is still much longer than £ d rag, so that 
the dust responds coherently to the eddy. Below this 
scale, we assume that the turbulence is subsonic and in- 
compressible so that it adopts a Kolmogorov spectrum. 
This assumption is conservative in the sense that contin- 
uing the I 1 / 2 scaling of supersonic turbulence to smaller 
scales would result in smaller grain velocities and thus a 
higher probability of dust sticking. For the Kolmogorov 
cascade, the velocities of the turbulent eddies scale like 
"oddy ~ c s (Vsonic) 1/3 , where l sonic is the length-scale of 
eddies with turnover times of t SO nic- The velocity of ed- 
dies with I < / son i c Or t < t S onic IS 



Veddy 



eddy 



1/3 



t-eddy 
^sonic 
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(14) 



The most important eddies are those with t e ddy ~ tdrag- 
For teddy <C ^drag, the dust does not respond to the eddy. 
In the opposite extreme, teddy tdrag, ah local dust 
particles will attain a net coherent motion. Hence the 
local random velocity of dust particles is set by the scale 
on which t od d y ~ Wag, i.e., 
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Numerical experiments indicate that dust grains that ap- 
proach one another at < 10 m/s will stick to each other 
(Poppe, Blum, & Henning 2000), while at much higher 
velocities ~ 1 kms -1 they are destroyed on impact. 

The thickness of the accretion disk in AGN is uncer- 
tain. Models predict h/R ~ 0.1 — 10 -3 on ~ parsec scales 
(e.g., Thompson et al. 2005), although this depends on 
the uncertain viscosity in the disk. Taking h/R ~ 10~ 2 , 
we find that the velocity dispersion for /im-sized dust 
grains is ~ 20 m s _1 . Due to the a 1 / 2 scaling for the 
random velocity from equation (|15j) , coagulation is more 
likely for smaller grains, while fragmentation becomes 
more likely for larger grains. In the previous section we 
showed that irradiation can extend the photosphere of 
the disk up to a height of ~ R. However, most of the 
mass of the disk is still contained in a thin layer with 
h/R <C 1 (see § 3). Thus the above argument about 
coagulation remains valid for most of the mass. 

We can estimate the timescale for dust coagulation by 
considering the fraction of dust grains which approach 
each other with a low enough velocity such that they 
will stick (e.g., Chakrabarti & McKee 2005). For this 
population, the coagulation timescale is 
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where / is the fraction of the dust particle population 
with v ~ 10ms -1 . The dust coagulation timescale is 
roughly the dynamical time for 1/im grains and much 
shorter for submicron-sized grains. Hence it is likely 
that a significant fraction of the dust grains in AGN 
disks will coagulate to form larger grains. As a check 
on our reasoning, equation (fTCl) can be scaled to condi- 
tions appropriate to proto-stellar disks. For M ~ 1M©, 
R ~ 200 A.U., and Q ~ 10 (Andrews & William 2007), 
we find tcoagAdyn ~ 1, similar to our result for AGN 
disks. Submm observations which probe dust emission 
from the outer parts of protostellar disks indeed find ev- 
idence for significant grain growth, consistent with the 
inference from equation (JTHJ) (Drainc 2006; Andrews & 
Williams 2007). 

In the ISM of the Milky Way, the timescale for dust to 
be destroyed by supernova shocks is significantly shorter 
than the timescale on which dust is injected by star for- 
mation (e.g., Jones et al. 1994). Thus grain growth in the 
interstellar medium appears crucial for establishing the 
observed properties of Galactic dust (e.g., the size distri- 
bution) . We suspect that a similar balance between grain 
growth (argued for above) and grain destruction will set 
the grain size properties in AGN disks. In the vicinity 
of AGN the intense radiation field, thermal sputtering 
by hot gas, and shock processing of grains all likely con- 
tribute to grain destruction, and preferentially destroy 
smaller dust grains (e.g., Laor & Draine 1993). Accu- 
rately estimating the grain destruction time is clearly 
difficult given the poor understanding of the physical 
conditions in parsec-scale AGN disks. As a plausible or- 
der of magnitude estimate we note that if the AGN disk 
is supersonically turbulent on scales where self-gravity is 
important (as in the Milky Way), dust will typically en- 
counter a supersonic shock once every dynamical time. 
During such encounters, dust grain are either broken up 
or have their outer parts vaporized. Freed metals quickly 
recondense (eq. [ID] ) and small dust particles coagulate 
(eq. [H]). Equation (fl"6|) implies that coagulation de- 
pends on dust grain size, with t coag ~ t^ yn for a ~ 1/im 
on parsec scales. If the dust destruction time is indeed 
~ idyn, this suggests that the AGN environment tends to 
favor the production of /im-sized dust grains, although 
this is clearly uncertain at the factor of a few level. 

If a significant fraction of the dust mass is contained 
in grains with sizes of ~ 1/xm the FIR opacity will be 
enhanced, increasing the importance of the irradiating X- 
ray and UV flux. For a spherical dust grain with radius 
a, the opacity is maximized when ira ~ A, where a is the 
size of the dust grain and A is the wavelength of interest. 
At T » 500 K, A « 6 fan and therefore the dust opacity is 
maximized for grains with a « 2/im. When the opacity 
is maximized, the cross section of the dust grain is close 
to geometric, so that the opacity is given by 

™ 2 fd oc . -i -i f fd \ 2 -i 

K ^i73w 25ai ^ 3 U^r mg ' ( r) 

where fd « 1% (for solar metallicity) is the dust mass to 
gas ratio, and as discussed above, increases roughly 
linearly with Z/Zq. Note that we have assumed a stan- 
dard dust grain density is p — 3gcm~ 3 . Large dust 
grains are likely full of holes and voids, lowering the den- 
sity and raising the opacity. Equation (|17[) implies a 
FIR opacity up to a factor of few-5 times larger than in 



the Milky Way (independent of any enhancement to the 
opacity due to supersolar metallicity). 

4. DISCUSSION AND CONCLUSIONS 

In the presence of a sufficient irradiating X-ray flux, 
the outer atmosphere of an initially thin AGN disk puffs 
up into a geometrically thick configuration that provides 
significant obscuration along most lines of sight. We esti- 
mate that for Lx > ix,crit ~ 0.01 — 0.1 LEdd X-ray heat- 
ing is dynamically important for the structure of AGN 
disks on parsec scales; the upper end of this range is ap- 
propriate for normal ISM dust, while the lower end of the 
range is appropriate if, as we have argued, the FIR opac- 
ity in AGN disks is enhanced through a combination of 
supersolar metallicities and grain growth (§4). The lat- 
ter may occur in AGN disks on parsec scales in direct 
analogy with massive protostellar disks. 

For Lx < ix,crit, irradiation has little dynamical effect 
on the properties of AGN disks, while for Lx > ix.crit) 
X-ray heating likely drives a significant outflow from the 
surface of the underlying accretion disk. The transition 
between these two different regimes is very abrupt with 
h/R increasing by a factor of ~ 100 with a factor of 
w 2.5 increase in X-ray flux (Fig. [3|). Our results thus 
predict a marked increase in the presence of small-scale 
nuclear obscuration for the most X-ray luminous AGN. 
This obscuring material should preferentially lie at dis- 
tances of ~ 0.1 — 3 pc from the central AGN because 

1 /2 

outside ~ 3Af 7 ' pc the bulge mass dominates the BH 
mass and so the irradiating flux is less likely to be locally 
super-Eddington. Grain growth, if it occurs, also occurs 
preferentially at radii of ~ 0.1 — 1 pc. Using the esti- 
mated mass loss rate in a continuum radiation pressure 
driven wind, we find characteristic obscuring columns of 
~ 10 22 " 23 cm" 2 for luminous AGN, in reasonable agree- 
ment with observations (e.g., Risaliti, Maiolino & Salvati 
1999). Without a more sophisticated model, however, 
it is difficult to quantitatively determine the solid angle 
subtended by this material, although we expect it to be 
~ 7r, implying a significant fraction of Type 2 AGN even 
at high luminosities. This is consistent with, e.g., the im- 
plied covering fraction of dusty material of ~ 10 — 50% 
from the IR fluxes of unobscured Type 1 quasars. 

Luminous Seyferts and quasars have X-ray Eddington 
ratios in the range required for X-ray irradiation to be 
dynamically important. This is particularly true if the 
FIR opacity is enhanced in AGN disks. In a sample of 
35 reverberation mapped AGN, Peterson et al. (2004) 
found that a significant fraction had LboI <^ O.lLEdd- 
Taking into account the X-ray emission up to ~ 100 keV, 
Seyferts typically have Lx/^Boi ~ 0.1 — 0.3 (e.g., Mar- 
coni et al. 2004), and so we expect quite a few systems 
will have X-ray fluxes that are super-Eddington with re- 
spect to the dust in a parsec-scale disk. At high red- 
shift, bolometric Eddington ratios for luminous AGN are 
w 1/3— 1 (e.g., Kollmeier et al. 2006), although quasars 
emit a smaller fraction of their flux in the X-ray than do 
Seyferts (e.g., Elvis et al. 1994; Marconi et al. 2004) and 
so their X-ray Eddington ratios are similar. 

We now discuss our model in the context of observa- 
tions of few specific systems. INTEGRAL observations 
directly probe the high energy emission from AGN that 
dominates the Compton heating of surrounding mate- 
rial. To focus on a concrete example, the Seyfert 2 
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galaxy NGC 4388 has an X-ray luminosity of Lx ~ 
4 x 10 43 ergss _1 between 20 — 100 keV and is Compton 
thin with N H w 3 x 10 23 cm" 2 (Beckman et al. 2006). 
NGC 4388 has a central BH mass of M BH « 6 x 10 6 M 
(Woo & Urry 2002), which implies L x /L Edd « 0.04. 
This X-ray Eddington ratio and the observed obscuring 
column are consistent with our results on the effect of 
radiation pressure on parsec-scale disks in AGN. 

To consider another concrete example, the prototyp- 
ical Seyfert 2 NGC 1068 is believed to be radiating 
at roughly a third of Eddington (Ho 2002; Pounds & 
Vaughan 2006). The inferred X-ray Eddington ratio for 
1068 is L x /L Edd « 0.01 (e.g., Panessa et al. 2006), 
but this is extremely uncertain because 1068 is Comp- 
ton thick. For reasonable intrinsic X-ray Eddington ra- 
tios, the X-ray flux of NGC 1068 is quite likely to be 
dynamically important for the surrounding material. In- 
frared interferometry of NGC 1068 at 8.7/zm by Jaffe et 
al. (2004) reveals the presence of a hot optically thick 
dust component (T > 800 K) 0.7 pc from the black hole 
with a scale height of h/R > 0.6. The spatial scale and 
geometric structure of this obscuring material are con- 
sistent with the predictions of our model. However, the 
obscuring column in 1068 exceeds our nominal prediction 
of ~ 10 22 — 10 23 cm~ 2 . In our models, larger columns 
would be observed at low inclination when the line of 
sight passes through the denser parts of the disks's at- 
mosphere (see Fig. 2}. 

Although radiation pressure is likely to be dynamically 
significant in shaping the observed properties of obscur- 
ing material in luminous AGN, it is equally clear that 
there are a number of systems in which this is not the 
case. In particular, low luminosity AGN such as NGC 
4258 (Greenhill et al. 1995) and NGC 3227 (Davies et 
al. 2006) are far too sub-Eddington for irradiation to be 
important. In such systems, other modes of obscuration 
are likely important, perhaps including a warped disk, 
obscuration by the host galaxy (e.g. MCG-6-30-15; Bal- 
lantyne, Weingartner, & Murray 2003), an accretion disk 
outflow, or inflowing molecular clouds. Our prediction 
is that luminous AGN should preferentially have a very 
compact source of nuclear obscuration that is intrinsi- 
cally geometrically thick (in contrast to a warped disk). 
Continued progress on IR interferometry should allow 
better observational probes of the spatial scale of nuclear 
obscuration in AGN. Continuum reverberation mapping 
between the IR and optical in a sample of Seyfert Is al- 
ready provides evidence that the obscuring material ex- 
tends down to the sublimation radius in many systems 
(e.g., Suganuma et al. 2006). 
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Our calculations have shown that geometrically thin 
disks on parsec scales in luminous AGN are likely to be 
substantially modified by X-ray irradiation. Once the 
surface layers of the disk start to become geometrically 
thick, UV heating will become dynamically important as 
well. Since the UV luminosity is larger than the X-ray 
luminosity, UV irradiation will ultimately have a signif- 
icant effect on the structure of the geometrically thick 
obscuring material (as in Pier & Krolik 1992), although 
the X-rays will likely, as we have argued, play the key role 
in initially thickening this material. A proper treatment 
of the structure of the irradiated parsec-scale disk thus 
ultimately demands a self-consistent multi-dimensional 
calculation that is beyond the scope of this work (see 
Krolik 2006 for axisymmetric hydrostatic calculations of 
radiation pressure supported tori). In light of these un- 
certainties, the mass outflow rate (eq. [5]) and obscuring 
column (eq. [9]) estimated in §3 can only provide an 
order of magnitude guide to the mass and column sup- 
ported by irradiation. Another over-simplification in our 
calculation is that we have modeled the disk as a uniform 
medium, which is probably not a good approximation. In 
addition to clumping due to self-gravity, radiation sup- 
ported atmospheres and winds are subject to instabili- 
ties that amplify inhomogeneities (e.g., Raylcigh- Taylor 
and photon bubble instabilities). This may increase the 
typical obscuring column at fixed M relative to that es- 
timated in equation Finally, we note that the mass 
outflow rate driven by radiation pressure (eq. [8]) may be 
substantial relative to the accretion rate onto the central 
black hole. This supports the hypothesis (e.g., Murray, 
Quataert, & Thompson 2005) that feedback from radia- 
tion pressure can help regulate the growth of supermas- 
sive black holes in galactic nuclei. 
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